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Chiral fluctuations in triangular antiferromagnets (TAFs) at T <^ Tm are studied theoretically. 
The case of a ferromagnetic interaction along c axis (which is directed perpendicular to the plane of 
the lattice), Dzyaloshinskii-Moriya interaction D||c and a weak magnetic field H||c is considered in 
detail. Previously, this model has been proposed to describe quantum TAP CsCuCls. Expressions 
for dynamical chirality (DC) are derived within the linear spin-wave approximation. In contrast 
to non-frustrated antiferromagnets, DC is found to be nonzero even at D, H — in a one-domain 
sample. We argue that this unusual behavior stems from the fact that a ground state of XY and 
Heisenberg TAPs is characterized by an axial vector along which DC is directed. 
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I. INTRODUCTION 

Spin chirality has attracted a lot of attention during last two decades (see Refs.-'^' ^i^i'^i^i^i'^ and references therein). 
The main subjects of investigation were helical magnets and triangular antiferromagnets (TAFs). A new (chiral) 
universahty class of phase transitions has been proposed for these systems by H. J. Kawamura (see Refi for a 
review). In particular, he observed that the dimension of their order parameter differs from those of antiferromagnets 
(AFs) on bipartite lattices. Let us discuss this point in detail for TAF. 

Since an order parameter space is defined as a topological space isomorphic to the set of ordered states)^ we consider 
the properties of the ground state. Due to frustration it has 120° spin structure. The ground state in this case can 
be characterized by chiralities of elementary triangles. The spin chirality is defined asi 

Cij — [Si X Sj], (1) 

where Si is a spin on cite i. The spin chirality of each elementary triangle is defined as a sum of chiralities of its spins 
when bypassing the triangle clockwise. The following two variants are possible: angles between all neighboring spins 
are a) 120° and b) 240° (see Figs. a) and b), respectively, for XY TAF). Chiralities of such triangles are directed 
oppositely. In Heisenberg TAFs spins not necessarily lie in the plane of the lattice. In this case the chirality of a 
triangle can have any direction but the spin configuration with the opposite chirality also exists. By convention we 
assume that triangle a) has positive chirality and triangle b) has the negative one. 

It is easy to show that there are two inequivalent 120° spin structures shown in Fig.^c) and d) which differ one 
from another by chiralities of each triangle. The spin order in configuration c) is determined by^ 

Sr = acos(koR) + bsin(koR), (2) 

where ~ ~ S^, (ab) — and kg is an antiferromagnetic vector which could be written in three forms (see 
Fig.Hl): 

C=2.(0.|0), ..". = 2,(-l,-i,„), .™=.,(-J,,-i,„), ,3) 

where the chemical lattice constant is taken equal to unity. As is seen from Fig. |2 they are equivalent up to the 
vectors of reciprocal chemical lattice. Spin configuration shown in Fig. ^ d) is described by Eq. jSJ with — kp put 
instead of kp. Subsequently, configurations c) and d) are referred to as (+ko) and (— ko), respectively. 

It is clear that XY TAF has two inequivalent ground states which differ by chiralities of each elementary triangle. 
This degeneracy of the ground state of XY TAF is called chiral degeneracy. As a result the symmetry of the order 
parameter is Z2XS1 instead of Si for coUinear XY AFs, where Z2 is the two-element group. In Heisenberg TAF the 
ground states (+ko) and (— ko) can be mutually transformed by global spin rotation. Nevertheless, the symmetry of 
the order parameter differs from that of non-frustrated AFs, as well. This is clear from the fact that the ground state 
in collinear AF is invariant under global spin rotation around magnetization of a sublattice whereas there is no such 
invariance in Heisenberg TAF. In Refi it is demonstrated that the symmetry of the order parameter in Heisenberg 
TAF is S0(3) instead of expected 83. 

Critical indexes of the chiral universality class calculated by Monte-Carlo simulations, e- and 1/n- expansions are 
found to be different from those of Heisenberg magnets, XY ferro- and XY antiferromagnets.^ The largest deviation 
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is in the index of specific heat: a = 0.24 ± 0.08 (n = 3) and a = 0.34 ± 0.06 in = 2) for the chiral universality class 
and a « —0.12 and a ~ —0.02, respectively, for usual magnets. These findings of the theory have been confirmed in 
numerous experiments^ 

It should be noted an important role of polarized neutron scattering in investigation of the spin chirality. Nowadays, 
it gives the only opportunity to verify one of the most intriguing results by Kawamura that the spin chirality is 
a new critical variable with fluctuations characterized by new critical indexes (3c and jc- The former describes the 
temperature behavior of the average C and the latter characterizes its fluctuations near the transition temperature. 
Index f3c can be investigated by elastic polarized neutron scattering.^' Experimental observation of 7c is a more 
complicated task. As C is a composite two-spin operator, the corresponding fluctuations are described by four-spin 
susceptibility and their observation is hardly possible. At the same time, as is shown in Refsi^, in presence of an 
axial vector interaction (e.g., Zeeman) the projection of C on it can be studied by polarized neutrons via the chiral 
term appearing in the cross section of inelastic scattering. 

Experimentally, the chiral critical indexes in XY TAFs have been determined recently for two members of the 
family of hexagonal AFs ABX3 {A — Rb, Cs; B = Mn, Fe, Co, Ni, Cu; X = CI, Br, I), where magnetic atoms form 
triangular lattice: CsMnBrg {S = 5/2)3 ^^^^ CsNiCla (S = 1)."^ The values of /3c = 0.44(2) and 7c = 0.84(7) obtained 
for CsMnBra and the sum /3c + 7c = 1.24(7) obtained for CsNiCla are in a good agreement with the theory. 

Thus, chiral fluctuations in TAFs near transition temperature have been extensively investigated recently and a 
certain success has been achieved. However chiral fluctuations in TAFs at temperatures much lower than the transition 
one have not been studied yet. One would expect that they have unusual properties in this case as well. 

In the present paper we consider chiral fluctuations in a quantum TAF CsCuCla at T <C TV. The spin- wave 
spectrum of this compound has been studied, e.g., in Ref^. However, the chiral susceptibility has not been addressed 
yet. Such an analysis would be important for future investigation of CsCuCla and other TAFs by polarized neutrons. 

CsCuCla is a unique member of the family of hexagonal AFs ABX3. Along c direction the copper atoms form 
ferromagnetic chains whilst in ab plane they interact antiferromagnetically and constitute a triangular lattice. The in- 
plane exchange is six times smaller than that along the chains. In contrast to CsMnBra and CsNiCls studied before, 
CsCuCls is characterized by Dzyaloshinskii-Moriya (DM) interaction along the chains. This compound exhibits a 
phase transition at Tjy = 10.7 K to the ordered state wherein spins lie in ab plane with 120° structure. They form 
a helical structure along the chains with a pitch 5.1° that is explained by a competition between ferromagnetic 
interaction within the chains and DM interaction.^^ We will assume also that a small magnetic field H||c is applied 
{H <C Hs, where Hs is a saturation field at which all spins become parallel). 

In this paper expressions for the dynamical chirality (or chiral vector) are derived within the linear spin-wave 
approximation. We find that in contrast to the non-frustrated AFs it is finite aX D,H = 0. As is well known^Sii^ the 
chiral vector is nonzero only in the presence of an axial vector (e.g., magnetization of the sample, DM interaction, 
etc.) along which it is directed. We demonstrate that a ground state of XY and Heisenberg TAFs is characterized by 
the axial vector. This is the consequence of existence of ground states (-l-ko) and (— ko). We show that these states 
are characterized by axial vectors [a x b] and — [a x b] , respectively. Hence, the chiral vector is perpendicular to the 
plane in which the spins lie and is directed oppositely for (-fkp) and (— kp) configurations. Therefore, the dynamical 
chirality can be observed in a sample with a different population of domains. 

The paper is organized as follows. Main principles of experimental investigation of chiral fluctuations by polarized 
neutron are considered in Sec.^] We present there expressions for the cross section, its chiral part and the dynamical 
chirality. We discuss general transformations of the Hamiltonian describing CsCuCla in Sec. IIIII Green's functions 
properties and spin-wave spectrum of this compound are studied in Sec. IIVI General properties of the chiral fluc- 
tuations in TAFs at T ^ T/v are discussed in Sec. IV Al Expressions for chiral vector in CsCuCls are derived in 
Sec. IV Bl within the linear spin-wave approximation. The range of validity of this approximation is also discussed 
there. Conclusions are made in Sec. IVII 



We consider in this section in some detail the main principles of experimental investigation of chiral fluctuations 
by polarized neutron. For the neutron scattering amplitude we have the well-known expressions 



II. INVESTIGATION OF CHIRAL FLUCTUATIONS BY POLARIZED NEUTRON 



Nq + Mq(T, 



(4) 




(5) 



(6) 
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where A'' is the total number of unit cells, &„ are nuclear scattering lengths, r = 5.39 x 10 cm, Fm{Q) are 
magnetic form- factors and Q = Q/Q, where Q is the momentum transfer. Four contributions to the inelastic 
scattering cross section can be represented using generalized retarded susceptibilities which have the formi^^ {A, B)aj = 
{A, B)'^ +i{A, -B)", where the first and the second terms are dispersive and absorptive parts, respectively. As a result 
for the cross section of inelastic scattering we have the following expressionii^ 

— - ^ [(A^-Q, A^q)" + (M-Q, Mq):: + zPo([M_Q X Mq]):: + Po(7V-q, Mq + M_Q, TVq):^] , 



dVLduj TT ki 1 — cxp(— w/T) 

. (7) 

where Pq is the initial neutron polarization. The first two terms in Eq. Q are nuclear and conventional magnetic 
parts of the cross section, respectively. The third and the fourth terms determine polarization-dependent part of the 
cross section that corresponds to chiral magnetic scattering and nuclear-magnetic interference, respectively^ As is 
seen from Eqs. © and |Q, one can eliminate contribution of the interference by directing polarization Pq along the 
impulse transfer Q. In this case the chiral part of the scattering is obtained experimentally by subtracting results of 
the cross section measurements with polarization Pq and — Pq. The chiral term in Eq. Q is proportional to imaginary 
part of dynamical chirality C(cij, Q) which is defined from the antisymmetric part of non-diagonal components of the 
generalized susceptibilityt^ 







where a, (3 = x^y, z and (...) denotes thermal average. Multiplying this expression by icajj-y we obtain the dynamical 
chirality: 



1 1"°° 

C{oj, Q) = - 2 '^^e^"* ([S-qW X Sq(0)] + [Sq(0) x S_Q(t)]) . 



(9) 



As a result we have for the chiral part of the cross section from Eqs. © and 



dflduj TT ki 1 — cxp(— cj/T) 



2r^F„^,(Q)(PoQ)(QImC(c.,Q)). (10) 



Then, the main object of discussion of the present paper is the chiral vector C{uj, Q). 
For the elastic scattering cross section we havet^ 

^ = (iV_Q)(iVg) + (M_q)(Mq) +zPo[(M_q) X (Mq)] +Po((A^_q)(Mq) + (M_g)(iVQ)). (11) 

One can express the chiral term in Eq. Hll|) via the staggered chiralityii 

C(Q) =i[Sq X S_q]. (12) 

We have from Eqs. ^ and 1)1111 that the chiral term in the elastic scattering cross section is proportional to C^(Q)Po, 
where 

C^(Q) - ( ([a X b] • q) ) 5][Z\(Q - ko + r) - A{Q + ko + r)], (13) 

T 

where r is the reciprocal lattice vector, A{0) = TV, Z\(Q 0) = and Q = Q/Q- Notice that the dynamical chirality 
given by Eq. is a natural generalization of the staggered one given by Eq. 1)12(1 . It is seen from Eq. ()13|l that 
C^(Q) is expressed via the average static chirality which can be investigated near Bragg reflections. Since Eq. H13|) 
is an odd function of ko, the polarization-dependent part of the cross section has opposite signs for states described 
by kg and — ko. However, if the energy does not depend on the sign of antiferromagnetic vector (as in the case of 
centrosymmetric crystals CsMnBrs, CsNiCls and helimagnetic holmium) the sample splits onto domains with ko and 
— ko. Therefore, the chiral scattering will be proportional to the difference of their population and can be observed 
only if this difference is nonzero. For CsMnBra it is the casei^ For helimagnetic holmium the difference in domain 
population has been achieved by cooling of the twisted crystal below the transition temperature^ 
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III. MODEL HAMILTONIAN 



The model Hamiltonian of CsCuCls has the iorm>— 
n = -2 Jo J2 (SmS„,+i + rj{SZ,SZ,+, + SlSl^,)) + 2Ji ^ S,„S,„ - ^ D • [S„ x S„+i] - gfigHj^SL, (14) 

in (ij)n in in 

where Si„ is operator of a spm at the i-th site in the n-th ab plane, z axis is taken to be parallel to c axis, Jq > and 
Ji > are values of exchange along c axis and within ah plane, respectively, < ^ 1 is a weak in-plane anisotropy, 
{ij) denote nearest-neighbor spins in ab plane. The third term in p4l) describes DM interaction with vector D being 
parallel to c axis. The fourth term is the Zeeman energy. Coefhcients in (|14|l have been estimated previously (see 
RefsiiiJA and references therein): Jq « 28 K, Ji « 4.9 K, D « 5 K and rj « 0.008. 

One can simplify the expression for Ti by applying a standard transformation of spin operatorst^i^ 



in 

•X' oiT, ^„ _L QV r.na ^r, ^ ' 



which implies a rotation of xy plane by a pitch q along z axis. The value of q is determined so as to eliminate the 
antisymmetric part of the Hamiltonian [Si„ x Sin+i]: 

D , , 

tang = (16) 

2Jo(l + ??) 

As a result the Hamiltonian has the form: 

= - ^ (2Jo(5,f^5,f^_,_i + "Sf^Sf^+i) + 2jQSinSin+ij + 2Ji ^ S-„S^„ - g^iBH^^ Sf^, (17) 

in {ij)n in 

where 



^0 = Joy (1 + ??)'+ «1.012Jo. (18) 

This way DM interaction has been reduced to a contribution to in-plane anisotropy. As is shown in Ref>ii, at = 
this anisotropy favors the coplanar ground state configuration with three sublattices and 120° spin structure. Spins 
are canted by magnetic field forming at iJ <C Hs an umbrella-like configuration. In the following derivations we 
assume this spin structure of the ground state and neglect the small difference between Jq and Jq. 

For consideration of spin waves it is convenient to use the following representations for the spin operators: 



= -/icos(koR™)Sf^' - sin(koR™)Sf" + x^^cos(koR„)5f^', 
Sf^ = -fe sin(koR ,„)5g;^ -I- cos(koR™)5f,;' + \/r~??sin(koR„)5,^', (19) 
S*?' — —\/l^^h?Sf'' — hS 



qxii _ S( , t <r,°-i,^ \ 

'^in ~ \ 2 V""i ^ "in 23 



QV"-^,,[I(a- -flt _^lAn\ (20) 
'-'in — 'A/ 2 v"*" "in 2S ' ' 



2 V"*" ^in 2S 

Qzll _ C _ „t „ . 
'^in ~ ^ "'in"«n: 

where /i > is a sinus of the canting angle [h — at H = 0), derived below, and ko is the antiferromagnetic 
vector. It should be stressed that the spin transformation (|19|) with antiferromagnetic vectors kp given by Eq. (0) 
describes configuration (-fko) (Fig.^c)). One has to use —kg instead of ko in p9|) to describe configuration (— ko) 
(Fig. Qd)). This difference has no effect on results of spin-wave spectrum calculations but, as is shown below, it is 
of great importance for chiral fluctuations. Unless otherwise specified, all formulas presented below are for (+ko) 
configuration. 

Note that representations (|19|) and (|20|l contain only one type of Bose operators. This makes intermediate calcu- 
lations and results more compact. Such an approach has been previously applied for other problemsiSiiS and proved 
to be more convenient than the traditional one dealing with several types of Bose operators ascribed to different 
sublattices. 



Substituting Eq. (|19|l into Eq. (|17() and taking the Fourier transform of spin operators: 

S» = \/5ESke^''''"' (21) 

k 

where N is the number of spins in the lattice and the sum is over the chemical Brillouin zone (BZ) (see Fig.O, we 
have for the Hamiltonian: 

7^ = -2Jo^cosfc,S;;S'ik + 6Ji^i^kS;^S'ik-gAiB^^\/^'S'^", (22) 

k k 

where 



f k = ( exp{ifcj,} + exp I i ( ^ ^Kj i + exp i j _ ifc^ ) [> ) . (23) 



Let us substitute Eq. (|20|l into Eq. ()22() . The classical ground state energy is given by the term without Bose operators: 

Eo = NS{gfiBhH - 2SJo - 3(1 - 3h'^)SJi). (24) 

Minimization of Eq. (|24|l with respect to h yields h = —giisH/Hs [Hs — ISS'Ji is the saturation field) in accordance 
with Ref^. In this case the terms linear in Bose operators cancel each other and the Hamiltonian has the form 
H. = Eq + X]l'=2^«' where Tii describes terms with products of i operators a and a^. Below we derive the spin 
chirality within the linear spin- wave approximation, i.e., we restrict ourself to bilinear part of the Hamiltonian which 
has a simple form: 



-^k 



El, = ASJoil - cosfc^) + 6SJi + 35Ji(l - 3h^)\k - eVsSJih^w, (25) 



where Ak — Ret'k and ^k = Imi^k- The range of validity of this approximation is discussed at the end of Sec. IV Bl 

IV. GREEN'S FUNCTIONS AND SPIN- WAVE SPECTRUM 

It is convenient to introduce two Green's functions: G'(w,k) = —{au,al.)ui and F{uj,'k) — — (ak,a-k)t^ (definition 
of {. ■ ■)uj is given in Eq. ©)• Spin-wave spectrum can be obtained from analysis of their denominator. There are 
two other Green's functions connected to the ones above: G{uj,'k) = — (aLk''^-k)w = G* {—uj, —k) and Ft((jj,k) — 
— (a^j^, aj^)i^ = F*{~u!, — k). As a result we have two sets of Dyson's equations. One of them has the following form: 



(26) 



G(w,k) = GW(cj,k) +G'(°H^,k)Bki^1'(t^,k), 
Ft(w,k) = G^°\uj,k)BkG{uj,k), 

where G'-"-' (w, k) — {u; — i?k + i5)^^ is the bare Green's function. Solving Eqs. H26() we have: 

G(a;,k) = ^^±%, G(c.,k) = ^^i±^, Ft(c.,k)=F(^,k) = --?t-, (27) 
A(a;,k) A(w,k) A(a;,k) 

where 

A{uj, k)^{Lu + i6f - (cj + i6){Ei, - £;_k) - ^k£^-k + Bl = {uj-ei, + iS){u; + e_k + iS). (28) 

Here 



ek = ^ [4S'Jo(l - cosfc^) -I- 6SJi{l - Ak)] [4S'Jo(l - cos k,) + 6SJi{l + (2 - Sh^)X]^}] - 6\/35 Ji/i^k (29) 
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determines the spin-wave spectrum which is known to have three branches. Eq. (|29|l describes all three branches at 
once because vector k in our consideration lies within the first BZ of the reciprocal chemical lattice which includes the 
first and the second BZ of the reciprocal magnetic lattice (RML) (see Fig.|21l. Consequently, the first branch of the 
spectrum is given by Eq. (|29|l with k within the first BZ of RML. The second and the third branches are described 
by Eq. H29() with k being within, respectively, dark colored and light colored areas of the second BZ of RML shown 
in Fig. 13 One can represent expressions for these branches in a traditional form writing 

(1) (2) (3) /gn^ 

where ek is given by Eq. H29|l and vector k lies within the first BZ of the RML. It can be easily verified using Eqs. ((SJ, 
(|23|l . (|25|) and H29|l that these expressions are in agreement with those of Refiii. 
In the vicinity of points k = 0, k = kg and k = — kg the spectrum has the form: 



et = < 



cj|'^«'fc,)'+(c^'"y (ki-ko)VA2 if |k-ko|«A;o, (31) 
c\^"h,)\ (c^-'-y (ki + ko)' if |k + kol « A:o, 

where c[|°^ = ^/2cf"^ = QS^/J^, cf^ = V2c^*^°^ = 3^35 Ji, k^ = {k^,ky,Q) and A = ISSJih. We see that the 

spectrum is linear near these points with different spin- wave velocities aX k <^ k^ and |k ± ko| <C /cq. Magnetic 
field leads to a gap A at k = kg. We neglect here a small gap related to DM interaction. In om description it is 
proportional to Jo — Jo (see Eq. H18(l 'l. 

V. CHIRAL FLUCTUATIONS 

We discuss in this section general properties of the chiral fluctuations in TAFs and derive the dynamical chirality 
for CsCuCls. 

A. General properties 

As is well known^Siii the chiral vector is nonzero only in the presence of an axial vector (e.g. magnetization of the 
sample, DM interaction etc.) and is directed along it. To illustrate this point we consider two examples: square 
antiferromagnet and ferromagnet. 

Chiral vector is zero in the first case. To show that Cy = Q {z axis is along magnetization of sublattices) we make 
a rotation of the coordinate space by an angle tt along x direction. The Hamiltonian does not change under this 
transformation and the ground state spin configuration is described by the same antiferromagnetic vector. Hence 
from Eq. Q we have Cy{oj, Q) = ~Cy{uj, Q) — 0. One can show in the same way that Cz = Cx = 0. 

In the case of a ferromagnet the chiral vector is directed along magnetization of the sample (one-domain case is 
considered for simplicity). It becomes clear if we direct z axis along magnetization m and make the rotation of the 
coordinate space by tt along x direction. The Hamiltonian does not change and the ground state wave function is 
transformed to that of state with magnetization of the opposite direction. Thus from ^ we have ci"^\uj,Q) = 
—Ci "^\lu,Q). Making the rotation along z direction it is easy to show that Cy = Cx — 0. Then, the chiral vector 
is directed along m. 

We demonstrate now that there is a different situation in XY and Heisenberg TAFs. Let us assume that the ground 
state has configuration (-l-ko) and z axis is directed perpendicular to the plane in which the spins lie. It is easy to see 
that after the rotation of the coordinate space by tt along x direction the Hamiltonian does not change and the ground 
state spin configuration becomes (— ko). As a result we have Ci'*'*'''^ (w, Q) — — Ci~''°^ (w, Q). The rotation along z 
direction g means that chiral vector in configuration (-l-ko) has the same value as and is directed 

oppositely to that in configuration (—kg). But it certainly does not mean that it is zero. Results of straightforward 
calculations for XY TAF presented below show that it is in fact nonzero. 

Thus we see that existence of inequivalent ground states described by Eq. Q with the same a and b and different 
antiferromagnetic vectors kg and —kg is the origin of this unusual properties of chiral fluctuations. As a result 
direction perpendicular to the plane in which spins lie is physically selected. It is clear that the staggered chirality at 
Q = ko can play the role of the axial vector characterizing the system. From Eqs. and (|12|l . C(ko) = i[a x b] for 
(-l-ko) and C(ko) = — ^[a x b] for (— ko) configurations. 
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It should be noted that in real substances there are domains with both spin configurations (+ko) and (— kp), so 
that we have to introduce concentrations 'T'(+ko) s-'^d 'T-(-ko) of these domains in the sample, with fi(+ko) + "-(-kp) = 1- 
If D,H = their contribution to the chiral vector has opposite signs and, similar to staggered chirality (|12|l . the 
dynamical one would be nonzero only when ?T.(+kn) "-(-ko)- 

It can be shown using symmetry with respect to time inversion and definition @ that the chiral vector in general 
has the following propertiesi^iii 

C(w,Q,H) = -C(l.;,-Q,-H), (32) 
ImC(tJ,Q,H) = ImC(-w,-Q,H) = -ImC(-tj,Q,-H). (33) 



B. Chiral vector in CsCuCla 

We proceed with calculation of the dynamical chirality in CsCuCla. As is shown in Sec. it is expressed via non- 
diagonal components of spin Green's function Q. In CsCuCls there is only one non-diagonal component, Xxy = —Xyx, 
and the chiral vector © is directed along z axis: 

C,iu;,Q)^^Xxyi^,Q)- (34) 

To calculate it within the linear spin- wave approximation one has to use Eqs. ^ and (|34() . perform sequentially 
transformations of spin operators ifT^ . ifT^ . and use expressions l|77|) for Green's functions. 

Let us consider firstly general expressions for the chiral vector. It is easy to show using Eq. (|21|l that after 
transformation (|15|l the spin operators have the form: 

1 7 

QX ^ / QXl , QXl \ : / CV CV \ 

-^Q — Q- Q+'' 2 ~ 
7 1 

/ QX' QXf \ I 

where 



S'a = -;^(^Q- - '^Q J + + S'^^), (35) 



Q+ = Q + qc, Q =Q qc. (36) 

Here q is given by Eq. and c is a unite vector directed along c axis. After performing transformation (|15|l . for 
the chiral vector we have: 

C.(w, Q) - i [xL(a;, Q-) - x'xxi^, Q+)] + ^ [x'xyi^, Q ) + x'xyi^, Q+)] , (37) 

where x^^(t^,Q) = (s'"q,S'q^ • Deriving Eq. ^ we also used xi::,(w, Q) = x'yyi^^M) and xi;y(t^, Q) = ~x'yxi^M)- 

Note that the first term here is nonzero only when q ^ (i.e. D ^ 0). In contrast to non- frustrated systems the 
second term in (|37|l and the chiral vector itself is nonzero a,t D,H ^ 0. After transformation H19|l for x'xxi'^i Q) and 
x'xyi'^iQ) appearing in H37|l we obtain: 

X'xxi^, Q) ^ \ [Xyyi^, Q + ko) + Xyyi^, Q " kp ) + {x'Li^ , Q + ko) + Xxxi^, Q " ^o)}] 

+ y [x'iyi^, Q + ko) - x^',(^, Q - ko)] (38) 

X'xyi^, Q) = ^ [Xyyi^, Q + ^o) ~ Xyyi^, Q " ko) + {x'U^, Q + ko) - Xxxi^, Q ^ ko)}] 

- ^ [Xxyi^, Q + ko) + Xxyi^, Q - ko)] , (39) 

where x'afsi^^ Q) = ^"^-'q' "^q"^ • derive an explicit expression for chiral vector we have to make transformation 
(|20|l and use Eq. (|27|) for Green's functions. To calculate the corresponding expressions for configuration (— ko) one 
has to replace kg by — kp in H38|l and (|39|l and note that bilinear part of the Hamiltonian for (— ko) has the form l|25|) 
with a formal replacement of hhy —h. 
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As a result of straightforward calculations, using n(_|_k^) = 1/2 + (n(+ko) — 'T-(-ko))/2 and ?T-(_ko) = 1/2 — (n(+ko) — 
n(_i(;o))/2 and noting that terms proportional to 'T.(+ko) ^ '^(-ko) cancel each other, we have for the chiral vector at 
Q < ko: 



1 



72 Ji 



36SJiUjh 



ZQSJiUjh 



A(w,Q_+ko) A(w,Q+-ko) 



(40) 



As is seen from Eq. (I40II . the chiral vector is zero a,t D,H = and Q <C fco. 
For IQ — kol <C fco one has: 



ISS-Ji 



ISS-Ji 



ISSJiUjh 



ISSJlLuh 



A(cj,Q_ - ko) 

"(-ko)] 



S 



A(cj,Q+-ko) 
ISSJi 



9SJiA{lo,Q- 
ISSJi 



95JiA(cj,Q+) 
+ ISSJiujh 



ISSJiLuh 



A{u;,Q--U) A(w,Q+-ko) 95 JiA(w, Q_ + ko) 



9SJiA{u;,Q+) 

(41) 

and gives a finite 



The first term in Eq. H41(l vanishes aX, D,H — 0. The second term is proportional to Ti(+ko) ~ 
contribution to chiral vector even ai D, H = provided that ?^(+ko) "-(-ko)- accordance with the results above 
at D,H = the domains of the type (+ko) and those of the type (—ko) give contributions with opposite signs. 
Finally, for |Q + ko| <C fco the chiral vector becomes: 



s 
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ISSJi 



ISSJi 



Q 



ISSJiLuh 



Q 



- ISSJiUjh 



A(w,Q_ 



*(+ko) 



A(w,Q+ + ko) 95JiA(w,Q- 



95JiA(cj,Q+ - ko) 



H-ko 



ISSJi 



ISSJi 



e-Q_ + ISSJiLuh e2_Q^ - ISSJiUjh 



A(w,Q_ 



A(cj,Qh 



9SJiA{uj,Q-) 95 JiA(w, Q+ - ko 



This expression has the same structure as that for |Q — ko| ^ /co discussed above. 

As a result we see that even at D, H = the chiral vector remains finite in the vicinity of points ko and 
long as the population of domains of the type (+ko) and (— ko) are different. 

Using expressions H4()ll - H42() and 



Im 



1 



A{uj,k) Ek + e-k 



[Siuj- 



e-k 



S{lli - ek)] 



(42) 



-ko, as 



(43) 



one can verify that relations H32|) and H33|) hold in our case. In particular, the terms in expressions for C and ImC 
with different /i-parity have different uj- and Q- parity. Indeed, according to H32|) terms in C which are odd functions 
of h are even functions of Q and vice versa: even-/i terms are odd functions of Q. Similarly, we determine from l|33() 
that even-h and odd-h terms in ImC are odd and even functions of lu, respectively. 

Using Eqs. H1U|) . H4U|I -H43 () it is easy to derive the chiral part of the cross section. Because of cumbersomeness of 
the resulting expressions we do not present them here. The only thing we note is that according to Eqs. H33II and p(J|l 
at a; ^ r even-h and odd-h terms in the cross section are odd and even functions of uj, respectively. 

Let us consider now the range of validity of final expressions for the chiral vector (|40|I - H42|I . To derive these 
expressions we used the linear spin-wave approximation restricting the initial Hamiltonian (|14|) after transformations 
(|15|l . (|19|l and H2U|) to its bilinear part in operators of creation and annihilation H25|l . This approximation is widely 
used because it is believed that 1/5-corrections stemming from high-order terms in Hamiltonian are small. It is the 
case, e.g., for square and cubic Heisenberg AFs.^^'^^'^^'^'''^"''^^ At the same time the spin- wave interaction is important 
in some cases. For example, in square and cubic AFs in the magnetic field the first 1/5-correction to the chiral vector 
become large at small T, uj and QM- We carried out the corresponding analysis for TAFs. We found that the first 
1/5-corrections to the spin- wave spectrum and to the chiral vector are small. Therefore, one can restrict ourselves by 
the results of linear spin- wave approximation H4U|) - (|42() in the range of validity of representation H2U|) . i.e., at T ^ T^. 
The detailed discussion of this question is out of the scope of this paper. 



VI. CONCLUSION 



In this paper we study chiral fluctuations in triangular antiferromagnets (TAFs) at T <C T/y. The case of TAF in 
magnetic field Hp (c is a unit vector directed perpendicular to the plain of the lattice) with Dzyaloshinskii-Moriya 
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interaction Dp is considered in detail. This model has been proposed previously to describe quantum TAF CsCuCls. 
Expressions for dynamical chirality (or chiral vector) given by Eq. © describing the chiral fluctuations are derived 
within the linear spin-wave approximation. We demonstrate that in contrast to non-frustrated systems it is nonzero 
even a.t D,H — 0. As is well known,^-^^ the chiral vector is nonzero only in the case of presence an axial vector (e.g., 
magnetization of the sample, DM interaction etc.) along which it is directed. We obtain that a ground state of XY 
and Heisenberg TAFs is characterized by such an axial vector. It is perpendicular to the plane in which spins lie 
and is directed oppositely for the states with different chiralities of elementary triangles (see Fig. Possibility of 
investigation of the chiral fluctuations by polarized neutron is also discussed. 



Acknowledgments 

I am grateful to S. V. Maleyev for fruitful discussions and for his censorious remarks. It is a pleasure to thank A. G. 
Yashenkin and D. V. Lebedev for their help in preparation of the manuscript for publication. This work was supported 
by Russian Science Support Foundation, RFBR (Grant Nos. SS-1671.2003.2, 03-02-17340 and 00-15-96814), Grant 
Goskontract 40.012.1.1.1149 and Russian Programs "Quantum Macrophysics" , "Collective and Quantum Effects in 
Condensed Matter" and " Neutron Research of Solids" . 



* Electronic address: [syromyat @t hd . pnpi. spb . ru| 

^ H. J. Kawamura, J. Phys.: Condens. Mat. 10, 4707 (1998). 
2 S. V. Maleyev, Phys. Usp. 45, 569 (2002). 

^ V. P. Plakhty, J. Kulda, D. Visser, E. V. Moskvin, and J. Wosnitza, Phys. Rev. Lett. 85, 3942 (2000). 

* V. P. Plakhty, S. V. Maleyev, J. Kulda, E. D. Visser, J. Wosnitza, E. V. Moskvin, T. Bruckel, and R. K. Kremer, Physica 
B 297, 60 (2001). 

^ H. J. Kawamura, Rhys. Rev. B 38, 4916 (1988). 
S. V. Maleyev, Phys. Rev. Lett. 75, 4682 (1995). 

D. H. Lee, J. D. Joannopoulos, J. W. Negele, and D. P. Landau, Phys. Rev. Lett. 52, 433 (1984). 

* N. D. Mermin, Rev. Mod. Phys. 51, 591 (1979). 

^ S. V. Maleyev, V. G. Bar'yakhtar, and R. A. Suris, Sov. Phys. Solid State 4, 2533 (1963). 
^° B. Blume, Phys. Rev. 130, 1670 (1963). 

T. Nikuni and H. Shiba, J. Phys. Soc. Jpn. 62, 3268 (1993). 

K. Adachi, N. Achiwa, and M. Mekata, J. Phys. Soc. Jpn. 49, 545 (1980). 
" S. V. Maleyev, Physica B 297, 67 (2001). 

S. W. Lowesey, Theory of Neutron Scattering by Condensed Matter (Oxford University Press, Oxford, 1987). 
^5 V. P. Plakhty, W. Schweika, T. Bruckel, J. Kulda, S. V. Gavrilov, L.-P. Regnault, and E. D. Visser, Phys. Rev. B 64, 

R100402 (2001). 

" H. Tanaka, U. Schotte, and K. Schotte, J. Phys. Soc. Jpn. 61, 1344 (1992). 
D. N. Aristov and S. V. Maleyev, Phys. Rev. B 62, R751 (2000). 

D. Petitgrand, S. V. Maleyev, P. Bourges, and A. S. Ivanov, Phys. Rev. B 59, 1079 (1999). 
A. V. Syromyatnikov and S. V. Maleyev, Phys. Rev. B 65, 012401 (2001). 
^° P. Kopietz, Phys. Rev. B 41, 9228 (1990). 

2^ G. E. Castilla and S. Chakravarty, Phys. Rev. B 43, 13687 (1991). 

C. M. CanaU, S. M. Girvin, and M. Wallin, Phys. Rev. B 45, 10131 (1992). 
J.-I. Igarashi, Phys. Rev. B 46, 10763 (1992). 

A. B. Harris, D. Kumar, B. L Halperin, and P. C. Hohenberg, Phys. Rev. B 3, 961 (1971). 
^5 S. Maleyev, Phys. Rev. Lett. 85, 3281 (2000). 



10 




c) d) 




FIG. 1; a) and b) 120° spin structures of an elementary triangle with positive and negative chirality, respectively; c) and 
d) Two inequivalent 120° spin structures of a triangular antiferromagnet which differ one from another by chiralities of each 
triangle. They are described by Eq. J^J with the same a and b and with antiferromagnetic vectors ko and — ko, respectively. 
These configurations are referred to in the text as (+ko) and (— ko), respectively. 




FIG. 2: The reciprocal chemical lattice (RCL) (open circles) and the reciprocal magnetic lattice (RML) (black dots). The 
inner (white) hexagon is the first Brillouin zone (BZ) of RML. The outer hexagon defines the edge of the second BZ of RML 
and, simultaneously, the first BZ of RCL. Antiferromagnetic vectors kQ^', kg^^ and kp^' are also presented. They are equivalent 
up to vectors of RCL. Also shown are two parts of the second BZ of RML which can be reduced to the first BZ by shifting on 
ko (dark colored region) and — ko (light colored region). 



